Magic wavelengths for lattice trapped Rubidium four-level active optical clock 
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After pumped from 5s 1 / 2 ground state to 6p x / 2 state, the population inversion between 6s 1 / 2 and 
5Pi/2,3/2 will t> e established for Rubidium four-level active optical clock. In this paper, we calculate AC 
Stark shift due to lattice trapping laser which dominates the frequency shift of clock transition in lattice 
trapped Rubidium four-level active optical clock. Several magic wavelengths are found that can form 
desired optical lattice trapping potential. By choosing a proper intensity and linewidth of trapping laser, 
the fractional frequency uncertainty of clock transition due to AC Stark shift of trapping laser, is estimated 
to be below 10" 18 . 

PACS numbers 06.30.Ft, 32.10.Dk, 32.60,+i. 



Time or frequency precise measurement is extremely 
important for fundamental science and technology. Op- 
tical clock is, undoubtedly, the best time or frequency 
measurement instrument implemented to date|[l|-[3t]. The 
most accurate optical clock is Aluminium single-ion clock 
which has achieved a fractional frequency uncertainty of 

8.6 x 10 13 [|JLL|. Although single-ion optical clock pos- 
sesses high accuracy, its stability is limited according to 
Allan variance, a y (r) oc 1///VJ1. Using neutral atoms 
trapped in optical lattice is one way to break this limitation 
for it can greatly increase the total number of atoms Q5J] . 
Recently, a Strontium optical lattice clock with 10 -16 frac- 
tional uncertainty has been reported[|2l Hfl, using so called 
"magic wavelength" trapping light to cancel the light shift 
of the two clock transition levels involved[5]. In such op- 
tical lattice system, neutral atoms are trapped in optical 
lattice and the trapped neutral atoms experience frequency 
shifts which arise from various of mechanics, i.e., trapping 
light shift, black-body shift, collision shift, first and second 
order Doppler shift, etc. In the Lamb-Dick regime [|6l lZD, 
collision shift and the first and second order Doppler shift 
are greatly suppressed since their movements are confined. 
Therefore, most of the time, the first two terms will domi- 
nate the frequency shift. 

However, both single-ion trapped optical clock and neu- 
tral atoms optical lattice clock mentioned above are work- 
ing in passive mode which is limited by the linewidth of 
interrogating laser. On the other hand, an active opti- 
cal clock[18r-llin will output laser directly whose frequency 
is determined by optical clock transition. Therefore, it 
can provide two distinguishing qualities, i.e., a super- 
narrow quantum-limited linewidth and extremely small 
cavity pulling shift[9]. Theoretically, a Strontium clock 
using thermal atomic beam implies a quantum-limited 
linewidth of 0.51 Hz and outshines state-of-art narrowest 

6.7 Hz of Hg single-ion clock and 1.5 Hz of Strontium op- 
tical lattice clock[9]. It promises to improve the stability 
of the best clocks by 2 orders of magnitude, and it can be 
extended to quenching laser and superradiant Iaser [|l2l - ll5ll . 



Lattice trapped Rubidium four-level active optical clock 
(Fig. [U is a new experimental scheme to realize active opti- 
cal clock. A population inversion exists between the clock 

v Pumping 



Lattice / 




6P1/2 



5Si/2 



6S1/2 

.P°ck|/VWWv 

5P 3 / 2 

5Pl/2 




Repumping 



FIG. 1 : Scheme of lattice trapped Rubidium four-level active op- 
tical clock, with 421.7 nm pumping light pumps neutral Rubid- 
ium from ground state 5s 1 / 2 to excited state 6p 1 / 2 and 780 nm, 
795 nm repumping light repump atoms back to the optical loop, a 
population inversion between 6s and 5p states formed with time 
evolution. And atoms are trapped in optical lattice clock at magic 
wavelength where light frequency shift of the clock transition 
arises from the trapping light is cancelled. 

transition 6si/ 2 — ^Pi/2,3/2 with 8.4% in 6s 1 / 2 and 3.3% in 
5^3/2 when pumping from 5s to 6pi/ 2 with a laser of 421.7 
nm wavelength and 318 mW/cm 2 intensity fl6h . The popu- 
lation inversion is formed after 10 -6 s, thus a steady photon 
number will be sustained with a bad cavity to form an ac- 
tive optical clock. Considering the Rayleigh scattering and 
Raman scattering in optical lattice, still a trap lifetime lager 
than 100 s is achievable 1 1 714 1 9n . such a long period is suffi- 
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cient for lattice trapped Rubidium four-level active optical 
clock. In this paper we derive several magic wavelengths 
for such a system. 

Neutral Rubidium can be cooled and trapped in optical 
lattice due to dipole gradient force fl20l - l22ll . The trapping 
light, meanwhile, causes light shift of energy levels which 
is affected by both the wavelengths and the intensities of 
trapping light. 

According to quantum perturbation theory, energy shift 
of arbitrary state can be expressed as follow, 



AE = ~(E) 
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where E is the electric field strength, a y ) is the polarizabil- 
ity of related states. Here, we consider the electric dipole 
interaction and neglect all higher-order perturbations. 

We consider the linearly polarized light to exclude the 
axial polarizability associated with light polarization and 
choose the orientation of the trapping light polarization 
along the magnetic z-axis, finally the polarizability expres- 
sion writes as ll23ll . 
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where orj'^ and Oi ± J AZ) represent scalar and tensor polar- 
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izabilities, and all variables and signs follow Ref. [23]. 

All reduced matrix elements (RMEs) and energy levels 
needed in our calculation are listed in Table U 

Part of RMEs are derived from the equation d 2 — 
3c 3 W(l + 2J)/(4j 3 ). Where d is reduced matrix ele- 
ment, W is the transition probability, J is the total angu- 
lar momentum quantum number. Above equation has used 
atomic units, where m e — h — e — = 1. 

To verify the validity of our calculation for rubidium 
atom we reproduce some calculations from Ref. H23II then 
we calculate the AC polarizabilities of 5pi/ 2 , 5p 3 / 2 and 
6si/ 2 using Eqs. ©. 

The polarizability of 6si/ 2 results from the summation 
of all the levels that have electric-dipole interactions with 
6si/ 2 . The transition frequencies np-Qs is readily to ob- 
tain from NIST website [25] and other related papers 02411 . 
the reduced electric-dipole moment values are from M. S. 
Safronova et al. ll26ll who calculated the reduced electric- 
dipole moments up to 8p-6s using a relativistic all-order 
method. Using these values to calculate the polarizability 
of 6s state is enough for a rough evaluation, and we will 
discuss the accuracy in the last section. 

Fig 12 and Fig |3] show the difference frequency shift be- 
tween 6S1/2 and 5pi/ 2 levels while the trapping light fre- 
quency varies from 550 nm to 650 nm, 1200 nm to 1600 
nm. The light shift of 5p 1 / 2 is dominated by the electric 
dipole interaction 5p 1 / 2 -8s 1 j 2 at wavelength 595 nm which 
is lower than the resonance wavelength. At wavelength 
612 nm, hp\/ 2 is greatly shifted by two electric dipole in- 
teraction 5pi/2-8si/2 and 5pi/ 2 -6d 3 / 2 . Two magic wave- 
lengths, 1342 nm and 1421 nm, are found between 1200 
nm and 1600 nm. 



TABLE I: Neutral Rubidium energy levels in wavelength (nm) 
and electric-dipole matrix elements in (a.u.). For 5pi/ 2 ,5ps/ 2 
states, transitions to 5s, 6s, 7s, 8s and to 4d, 5d, 6d are included, 
and for 6s 1/2 state, transitions to 5p, 6p, 7p, 8p are included. 



transition 


wavelength (nm) 


RMEs (a.u.) 


5pi/2 - 5s 1/2 


794.98 ° 


4.231° 




794.979 * 


4.221 c 


5pi/ 2 - 6s 1/2 


1323.88 ° 


4.146° 


5P1/2 - 7s l/2 


728.20 ° 


0.953° 


5pi/ 2 - 8s 1/2 


607.24 ° 


0.502° 


5P1/2 - 4d 3/2 


1475.65 ° 


8.051° 


5pi/ 2 - 5d 3 /2 


762.10 ° 


1.35 ° 


5pi/ 2 - 6d 3 / 2 


S OA CL 

620.80 


1.07 ° 


5p 3 / 2 - 5s 1/2 


780.24 


5.977 




780.241* 


5.956 c 


5p 3 / 2 - 6si/ 2 


1366.87° 


6.05° 


5 P3/2 - 7s l/2 


741 ° 


1.35° 


$P3/2 - 8s l/2 


OlO. 1 D 


u. / wo 


°P3/2 ^ u 3/2 


1529.26° 


3.63 ° 


blDo In — 4:d,K In 
u J^3/2 ^^5/2 


1529.37° 


10.9° 


$P3/2 ~ $ d 3/2 


776.16 ° 


0.67° 


5p 3/2 - 5d 5/2 


775.98 ° 


1.98° 


5p 3 /9 — 6do/9 


630.1 ° 


0.51 ° 


5P3/2 ~ ® d 5/2 


630.01 ° 


1.51° 


6s 1/2 - 5pi/ 2 


1323.88* 


4.119 c 


6s 1/2 - 6p 1/2 


2791.29* 


9.684 c 


6*1/2 - 7 Pl/2 


1298.28* 


0.999 c 


6s 1/2 - 8pi/ 2 


1030.67* 


0.393 c 


6s 1/2 - 5p 3/2 


1366.87* 


6.013 c 


6s 1/2 - 6p 3/2 


2732.18* 


13.592 c 


6s l/2 - 7 P3/2 


1292.39* 


1.54 c 


6s 1/2 - 8p 3/2 


1028.67* 


0.628 c 



a Ref.|23|] 

*Ref.|'24,25 

c Ref.|H] 



FigH] and Figj5] show the difference frequency shift be- 
tween 6si/ 2 and 5p 3 / 2 levels while the trapping light fre- 
quency varies from 550 nm to 650 nm, 1200 nm to 1600 
nm. With respect to state 5p 3 / 2 (\rrij\ = |), the wave- 
length 610 nm is lower than the resonance 5p 3 / 2 -8si/ 2 
wavelength, while wavelength 624 nm is between the res- 
onance 5ps/ 2 -8si/ 2 and resonance 5p 3 / 2 -6d 3 / 2 . For state 
5ps/ 2 (\rrij I = |), the wavelength 575 nm lies on the blue 
side of the resonance 5^ 3 / 2 -6d 3 /2,5/2. For optical clock 
transition 6si/ 2 -5ps/ 2 (\rrij\ = |), 1331 nm and 1461 nm 
are magic wavelengths; for 6s 1 / 2 -5p 3 / 2 (\rrij\ = |), two 
magic wavelengths 1336 nm and 1453 nm are found. 

Magic wavelengths in the range of 550 nm to 650 nm 
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FIG. 2: Figure shows, respectively, the frequency shift of the 
states 5p]_/ 2 and 6s x / 2 - $Pi/2 is plotted in dash dot line (red), 

6s 1/2 is plotted in solid line (blue), and the different frequency bp 3/2 (mj = ±±) is plotted in dash dot line (red), 5p 3/2 (r 



FIG. 4: Figure shows, respectively, the frequency shift of 
the states 5p 3/ / 2 (raj = ±2), 5P3/2( m j — =•=§) an d ® s i/2- 



shift between them is plotted in short dot line (black). The dif- 
ferent frequency shift goes to zero at magic wavelength, which is 
marked as cross symbol. This figure is plotted with the assump- 
tion that the intensity of trapping laser is 10 kW/cm 2 . 
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FIG. 3: Figure shows, respectively, the frequency shift of the 
states 5pi/2 and 6s 1 / 2 . §Pi/2 is plotted in dash dot line (red), 
6s]_/2 is plotted in solid line (blue). This figure is plotted with the 
assumption that the intensity of trapping laser is 10 kW/cm 2 . 



and 1200 nm to 1600 nm are listed in Tableffll The fre- 
quency uncertainties at each magic wavelengths are listed 
in TabkHIIlby properly adjusting trapping laser intensities, 
i.e., trapping depths. 

The reason why choosing four-level quantum system 
rather than three-level, optical lattice clock rather than ther- 
mal atomic beam are as follow. As for thermal atomic 
beam clock, residual Doppler shift will be extremely large, 
however optical lattice clock remove that effect otherwise 
since atoms are laser cooled and trapped in Lamb-Dick 
regime. On the other hand, three-level system is affected 
by light shift of pumping laser. Four-level quantum sys- 
tem, however, can avoid such disadvantage by choosing 
clock transition without involving ground state which is 
connected to pumping laser directly. 



=b|) is plotted in dot line (blue), and 6s 1 / 2 is plotted in solid line 
(black). Light shift of the states 5p 3 / 2 (raj = ±^), and 6s 1 / 2 in- 
tersect at the magic wavelength which is marked as circle symbol. 
Triangle symbol marks the cross point of light frequency shift of 
states 5p 3 / 2 (raj — ±§) and 6s 1 / 2 . This figure is plotted with the 
assumption that the intensity of trapping laser is 10 kW/cm 2 . 
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FIG. 5: Figure shows, respectively, the frequency shift of 
the states 5p 3 / 2 (raj = ±2), 5p3/2( m j = =•=§) an d ® s i/2- 
5p3/2( m j — ±3) is plotted in dash dot line (red), 5p 3 / 2 (raj = 
±|) is plotted in dot line (blue), and 6s 1 / 2 is plotted in solid line 
(black). Light shift of the states 5p 3 / 2 (raj = ±^), and 6s 1 / 2 in- 
tersect at the magic wavelength which is marked as circle symbol. 
Triangle symbol marks the cross point of light frequency shift of 



states 5p 3 / 2 (r 



i and 6si/2- This figure is plotted with the 



assumption that the intensity of trapping laser is 10 kW/cm 2 . 



In our calculation, we used the reduced matrix element 
data up to 8p-6s transition. Since the frequency shift of 65 
due to 8pi/2 and 8p 3 / 2 states contribute less than 0.081%, 
at all magic wavelengths, to the total polarizability of 65, 
therefore we believe that states above 8p will contribute 
fewer. As for the polarizability of 5pi/ 2 , states 5s 1 / 2 , 6si/ 2 
and 4d 3 / 2 will dominate, states 8si/ 2 and 6d 3 / 2 totally 
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TABLE II: The magic wavelengths for optical clock transitions 
6si/2~5pi j2 an d 6s 1//2 -5p3 /2 and lattice trapping depths to the in- 
tensity of lattice laser Av/I (kHz/kW/cm 2 ) refer to those magic 
wavelengths are listed in this table. 



6*1/2 - 5pi/ 2 
A (nm) Ai//J (kHz/kW/cm 2 ) 


6*1/2 - 5p 3 /2 
A (nm) Az//J (kHz/kW/cm 2 ) 


595 -74.5 


575 


-69.7 


612 -78.2 


610 


-77.7 


1342 1196.7 


624 


-81.0 


1421 -3707.2 


1331 


-4453.4 




1336 


-995.2 




1453 


-2878.1 




1461 


-2753.4 



have less than 0.639% effect. 

In conclusion, four magic wavelengths 595 nm, 612 nm, 
1342 nm and 1421 nm are available for transition 6si/ 2 - 
5pi/ 2 , seven magic wavelengths 575 nm, 610 nm, 624 nm, 
1331 nm, 1336 nm, 1453 nm and 1461 nm are candidates 
for transition 6si/ 2 -5£> 3 /2 (listed in TableHJ). To obtain a 
10 -18 fractional uncertainty or even less, one must well 
control the linewidth of lattice laser, i.e., below 10 kHz as 
listed in TableHTfl With respect to magic wavelengths be- 
tween 550 nm and 650 nm, a further discussion is needed 
since those wavelengths are going to above the ionization 
limit of 6S1/2 (741 nm). In addition, dynamic stark shift in- 
duced by black-body radiation will be discussed elsewhere. 

TABLE III: Adjusting the intensities / (W/cm 2 ) of trapping laser 
to provide sufficient trapping depths Av (kHz), for each magic 
wavelengths A (nm) the frequency uncertainties (mHz) due to lat- 
tice laser frequency fluctuation (assumed to be 10 kHz) are listed 
in this table. 

A (nm) I (W/cm 2 ) Au (kHz) uncertainty (mHz) 



595 


1E4 


-745 


4.3E-4 


612 


1E4 


-782 


4.6E-4 


1342 


1E3 


1196.7 


3.2E-1 


1421 


200 


-741.4 


9.0E-3 


575 


1E4 


-697 


3.9E-4 


610 


1E4 


-777 


4.6E-4 


624 


1E4 


-810 


4.8E-4 


1331 


200 


-890.7 


2.0E-1 


1336 


1E3 


-995.2 


4.2E-1 


1453 


250 


-863.4 


4.6E-3 


1461 


250 


-688.4 


3.8E-3 



Moreover, the pumping light of wavelength 421.7 nm 
and intensity 318 mW/cm 2 contributing frequency correc- 
tions to optical clock transitions of 6si/ 2 -5pi/ 2 , 6si/ 2 - 
5P3/2( mj =±%)> 6si/2-5p 3 /2( m .=±!) are 0.45 Hz, 0.65 Hz 
and 0.28 Hz each and related uncertainties will be no more 
than 1 mHz when an order of 10 -3 intensity fluctuation of 



pumping laser is considered. Although our results for Ru- 
bidium atom is specified, it is ready to be generalized to 
other alkaline atoms. 
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